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COMPUTATION OF LOW SPEED CAVITY NOISE 


Ching Y. Loh* 

Taitech, Inc. 

Beavercreek, Ohio 45430 


Abstract 

Cavity noise is an important noise source. It is 
generated primarily by wall-jet shear layer and 
vortex shedding and an aeroacoustic feedback 
loop. At low flow speed, the situation is further 
complicated by the presence of the Helmholtz 
resonance of the cavity. 

In this paper, the space-time conservation ele- 
ment and solution element (CE/SE) method is 
employed to numerically study the near-field 
noise of a cavity at low flow speeds. The nu- 
merical approach is of the MILES type (mono- 
tonically integrated large eddy simulation) [4- 
6], The computed frequencies and sound pres- 
sure levels (SPL) compare well with the exper- 
imental findings. 

1 Introduction 

Noise prediction and noise control are important for 
aerospace and automobile industries. The noises due to 
aeronautics at high speed, such as airfoil slat, landing 
gear and openings, etc. are some of the major noise 
sources in aeroacoustics. For ground transportation at 
low Mach number, the automobile industry is concerned 
with the noises produced by the door gaps, the side mir- 
ror and the aerial and how they affect the comfort in 
the car. All these noises may be categorized as airframe 
noises that are generated by the interaction between the 
vortex streets in the turbulent wake or between the vor- 
tices and the solid body edge(s). The phenomena are 
further complicated by a possible aeroacoustic feedback 
loop or a possible Helmholtz fluid resonance. 

Cavity noise is one of the most important airframe 
noises. When flow passes over a cavity or opening, due 
to vortex shedding at the upstream edge of the cavity, in- 
tensive tone noises may be generated. Strong tonal oscil- 
lations occur in a feedback loop between the two edges 
of the cavity opening. First, the vortices generated and 
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shed from the upstream edge of the cavity convect down- 
stream, impinge on the other edge, and produce acoustic 
waves. Then, as the acoustic waves propagate either in- 
side or outside of the cavity to the upstream edge, where 
receptivity of the wall jet shear layer is high, a new insta- 
bility wave (vortex) is stimulated and shed, a feedback 
cycle is thus completed. The tone noise generated in 
this way may be categorized as fluid-dynamic oscillation 
noise. Most high speed flow (supersonic, transonic, or 
high subsonic flows) noises are generated by the fluid dy- 
namic feedback oscillation mode. However, at low flow 
speed [1 - 3], depending on the geometry of the cavity, 
another type of tone noise due to fluid resonant oscilla- 
tion within the cavity may occur. The sound waves inside 
the cavity may be longitudinal or transverse depending 
on the aspect ratio of the cavity. For a deep cavity at 
low flow speed, both major tone noise generating mech- 
anisms may coexist. 

Numerical investigations of cavity noise based on di- 
rect numerical simulation (DNS) need very fine grids and 
are hence expensive. Large eddy simulations (LES) al- 
low a coarser grid since only large eddies are resolved 
and the smaller ones are modeled by a subgrid scale 
(SGS) model. 

In the present paper, a MILES (monotonically inte- 
grated large eddy simulation, [ 4 - 6]) type approach is 
used to investigate a typical low speed cavity noise prob- 
lem. In all the cases, the grid may not be fine enough for 
the Kolmogorov scale of turbulence but is fine enough for 
the large turbulence structures (large eddies or instability 
waves ) associated with the tone noise frequencies. The 
MILES approach appears somewhat similar to LES, but 
there is no explicit filtering since the cell-averaging pro- 
cess is already equivalent to some spatial filtering. Due 
to certain built-in numerical dissipation in a MILES type 
finite volume scheme, the SGS (subgrid scale) model 
may not be needed. The recent space-time conserva- 
tion element and solution element method (CE/SE) [7- 
8] is a MILES type finite volume method with generally 
less numerical dissipation and is adopted for the com- 
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putation. As demonstrated in the previous papers, the 
CE/SE scheme is well suited for aeroacoustics computa- 
tion [9,13], Because of the CE/SE non-reflecting bound- 
ary conditions (NRBC), which are based on the first prin- 
ciple of plane wave propagation [14], a smaller near field 
computational domain can be used in the present numer- 
ical simulation and helps to save both memory and CPU 
time. 

The 2-D unstructured-grid Navier-Stokes (N-S) 
CE/SE scheme used here is briefly discussed in Section 
2. Section 3 illustrates the noise problem of a rectangu- 
lar cavity with an overhang at low flow speed, including 
the initial and boundary conditions. The numerical re- 
sults are presented and compared to available experimen- 
tal findings [15] in Section 4. They are then compared to 
the experiment data. Concluding remarks are drawn in 
Section 5. 


Ua = p /{ 7 - 1 ) + p(u 2 + v 2 )/2. 

The flux vectors are further split into inviscid and viscous 
fluxes: 

F = F ; — F V ,G = G ; — G v , 

where the inviscid fluxes are the same as in the Euler 
equations: 

F n = U 2 , 

F i2 = (7 - 1 )Ua+ [(3 - 7 )U$ - (7 - 1 )Ul] /2U U 
F % 3 = U 2 U 3 /U U 

Fa = iU 2 Ua/Ua - (7 - 1 )U 2 [U 2 + U 2 ] /2U 2 , 

Gi i = U 3 , G i2 = U 2 U 3 /U U 
G i3 = (7 - 1 )Ua+ [(3 - 7 )Ul - (7 - F)U 2 \ / 2 U u 


2 The Numerical Scheme 

In the present 2-D MILES approach, the space- 
time conservation element and solution element, CE/SE 
method is chosen. The CE/SE scheme is a low order 
(second order accurate in time and space) scheme with 
low dissipation. A detailed description of the method 
can be found in [7- 8], In this Section, a brief sketch of 
the scheme is given. 


GiA = iUzUa/Ua - (7 - 1)173 [Ul + U 2 \ / 2 U 2 , 

and the viscous fluxes are: 

F v i = 0, F v 2 = p(2 u x - |V-V), 

F v 3 = p(v x + Uy), 

F v 4 = p[2 uu x + ( u y + v x )v — |(V*V)u+ 


2.1 Conservation form of the 2-D unsteady 
compressible Navier-Stokes equations 

In general, the CE/SE method systematically solves a 
set of integral equations derived directly from the physi- 
cal conservation laws and naturally captures shocks and 
other discontinuities in the flow. Both conservative vari- 
ables and their derivatives are computed simultaneously 
as unknowns. 

Consider a dimensionless conservation form of the un- 
steady 2-D Navier-Stokes equations of a perfect gas. Let 
p, u, v, p, and 7 be the density, streamwise velocity com- 
ponent, transverse velocity component, static pressure, 
and constant specific heat ratio, respectively. The 2-D 
Navier-Stokes equations then can be written in the fol- 
lowing vector form: 


U t + F x + G y — 0, ( 1 ) 


where x, y, and t are the streamwise and transverse co- 
ordinates and time, respectively. The conservative flow 
variable vector U and the flux vectors in the streamwise 
and the normal (y) directions, F and G, are given by: 


U = 


( U 1 

u 2 




G = 



7 d Ua u 2 + v 2 

2 

G„i = 0, G v 2 = p(v x + Uy), 

G v 3 = p(2 Vy - §V-V), 

GvA = p[2VVy + (Uy + v x )u — |(V-V)u + 

7 d Ua u 2 + u 2 

where u, v, u x , u y ,v x , v y are respectively the x — and y— 
flow velocity components and their derivatives, which 
can be written in terms of the conservative variables 
Ui, U 2 , U 3 and Ua , with Pr (= 0.72) being the Prandtl 
number, // the viscosity, the velocity divergence 

V*V = u x + Vy. 

By considering (x, y, t) as coordinates of a three- 
dimensional Euclidean space, E 3 , and using the Gauss 
divergence theorem, it follows that Eq. (1) is equivalent 
to the following integral conservation law: 

<f H m • dS = 0, m = 1,2,3, 4, (2) 

Jstv) 


with 


U\ = p, U 2 = pu , U 3 = pv, 


where S'(U) denotes the surface around a volume V in 

Pi and H m = (F m , G mi Um). 
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hexagon cylinder ADBECFA - A'D'B'E'C’F’A' 
wi1hte3QEs:0ADB-0’AD’B’, OBEC-O’B’E’C’, 
and OCFA - O’CFA’ 


Figure 1: CE/SE unstructured grid in space-time E$ 
space. 

2.2 Updating with a compact cell stencil 

The CE/SE scheme is constructed to take advantage of 
an unstructured triangle grid. The unstructured geome- 
try used with the CE/SE scheme is illustrated in Fig. 1. 
Here, A ABC is a typical triangular cell centered at O 
(or O' at the new time level) and D , E, F are the trian- 
gle centers of the neighboring cells. The flow variables 
at the previous time step are stored at the triangle cell 
centers. Centers D, E, F form a compact stencil to sup- 
port the marching of U at O to the next time level O' . 
Three quadrilateral cylinders (conservation elements) are 
formed by the edges that connect the vertices and the 
center of the triangle and its three neighbors. In the 
space-time E : >, space, Eq. (2) is applied to the hexagon 
cylinder ADBECF — A'D'B'E'C'F' of volume V that 
consists of these 3 quadrilateral cylinder CE s, see Fig 1. 

In the CE/SE scheme, the above flux conservation re- 
lation, Eq. (2), in space-time is the only mechanism that 
transfers information between node points. A conserva- 
tion element CE (here, quadrilateral cylinders) is the fi- 
nite volume to which Eq. (2) is applied. Discontinuities 
are allowed to occur in a conservation element. A so- 
lution element SE associated with a grid node (e.g., D, 
E, or F in Fig. 1) is here a set of interface planes in E$ 
that passes through this node (e.g. DAA'D ' , DBB'D', 
EBB'E', ECC'E' , etc.). Each surface S(CE) is made 
up of segments belonging to two neighboring CE’ s. 
Within a given solution element SE(j, n), where j, n are 
the node index and time step, respectively, the flow vari- 
ables are not only considered continuous but are also ap- 
proximated by the linear Taylor series expansions. The 


surface flux can then be calculated accurately and eas- 
ily by first evaluating the flux vectors at the geometrical 
center of the surface through such Taylor expansions. 

At time level n, the solution variables U, U x , and 
U y are given at the three nodes D, E, F in Fig. 1 and 
U, U x and U y at O' at the new time level n + 1 are 
to be computed. In principle, each of the 3 CE s pro- 
vides 4 scalar equations when Eq. (2) is applied to the el- 
ement. Hence, the 12 scalar equations needed for the 12 
scalar unknowns at O' are available. All the unknowns 
are computed based on these relations. No extrapolations 
(interpolations) across a stencil of cells are needed. 

More details about the unstructured CE/SE method 
can be found in [8]. 

2.3 Non-reflecting boundary condition 
(NRBC) for CE/SE 

In the CE/SE scheme, NRBCs are constructed based on 
the first principle of plane wave propagation and the hy- 
perbolicity of the equation system [14]. There are vari- 
ous variants of the NRBCs [9-14], the following are the 
ones employed in this paper. 

For a grid node ( j , n) lying at the outer border of the 
domain, where j is the grid node index number and n the 
time step number, the Type I NRBC requires that 

(U X )] = (Uy)] = 0, 

while U] is kept fixed at the initially given steady bound- 
ary value. At the outflow boundary, where there are sub- 
stantial gradients in the y direction, the Type II NRBC 
requires that 

(U x )] = 0, 

while U] and (Uy)] are now defined by simple extrap- 
olation from the nearest interior node j', i. e . , 

U] = U], (Uy)] = (Uy)],. 

As will be observed later, these NRBCs are robust 
enough to allow a near field computation without dis- 
turbing or distorting the flow and acoustic fields. 

3 Low-speed Tone Noise of a Cavity with 
Overhang 

For high speed flows ranging from high subsonic, 
transonic to supersonic, typical cavity noise examples 
are the airplane wheel well and weapon bay noises. 
Therefore, cavity noises are often studied at high flow 
speed when aeroacoustic feedback loop plays a dominant 
role. Rossiter [1] has a semi-empirical model to describe 
the convection of the vortices and hence the feedback 
tone frequencies for high subsonic flows. However, for 
ground transport vehicles at low subsonic Mach numbers 
(typically around or below M=0.1), or in the presence of 
a deep cavity, the noise due to flow-resonant phenomena 
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Figure 3: Trangulated grid around the overhang of the 
cavity. 


Figure 2: auto-door gap noise problem and the cavity 
dimensions. 


occurs as well. Howe [2] analyzes the cavity noise at 
low Mach numbers. Chatellier et al. [3] have detailed 
analysis and experimental data for a rectangular cavity. 

In the present paper, for low speed flows with Mach 
number as low as 0.08, a cavity with an overhang at the 
opening (Fig. 2) is investigated. This is a prototype of 
the automobile door gap noise problem [15]. 

Figure 2 illustrates the geometric configuration and 
the actual dimensions of the cavity noise problem. Due 
to the presence of the overhang tongue, the two os- 
cillation modes, namely, the fluid dynamic oscillation 
due to aeroacoustic feedback, and the Helmholtz fluid 
resonance in the cavity may emerge at the same time 
and compete with each other. The problem is given as 
a benchmark problem at the 3rd CAA (Computational 
Aeroacoustics) Workshop [15]. In Fig. 2, C is the length 
of the cavity opening. 2C is then chosen as the length 
scale. The speed of sound (343 m/s ) and density of 
the ambient flow are respectively the scales for velocity 
and density. In the current computation, there are about 
85300 triangular cells in the computational domain. The 
domain ranges —1 ^ x ^ 2.1 and 0 ^ y ^ 3.2. 
The triangles are actually obtained by dividing a rect- 
angular structured mesh cell into 4 pieces. The rect- 
angular cell keeps a uniform size of Ax’ = 0.00625 
and Ay = 0.00625 around the area of the opening and 
the interior of the cavity (Fig. 3), but grows larger near 
the outer boundaries. A nondimensional time step size 
At = 0.0025, and a weighted average (van Albada lim- 
iter) index a = 0 are employed. 


3.1 Initial and boundary conditions 

Initially, the flow of the entire domain is set at the ambi- 
ent conditions, i.e., 

Po = 1-0, m o = 0, ^o = 0, po = 1/7) 

All spatial derivatives are set to zero. 

At the inflow boundary, the flow variables are specified 
as the uniform ‘plug’ flow: 

Pi = 1.0, Ui = M, Vi =0, pi = I/ 7 , 

and their spatial derivatives are conveniently set to zero. 
The Mach numbers range from M = 0.0788 to M = 
0.147, corresponding respectively to flow speeds of 26.8 
m/s and 50 m/s. When a consistency test is conducted 
in §4.2.4, a prescribed profile is imposed on n-velocity at 
the inlet boundary. 

The no-slip wall boundary condition is applied to all 
the solid walls of the cavity body. Boundary condition 
at the top is the same as the inflow, which also acts as a 
Type I NRBC. At the outflow boundary. Type II NRBC 
is imposed. 

4 Numerical Results 

In this section, numerical results for the low speed cav- 
ity noises are presented and compared to their respective 
experimental data. As designed in the experiment, all the 
data were obtained at a location at the center of the inner 
left vertical wall of the cavity. 

As described earlier, for the overhung cavity consid- 
ered here, there are two major modes for the tone noises: 
the fluid dynamic tone noise due to the feedback loop 
around the cavity opening, and the fluid-resonant tone 
noises due to the Helmholtz resonance of the cavity. In 
addition, there could exist other tone noise modes of un- 
certain origins or broadband noises, making the phenom- 
ena rather complicated. For MILES approach, Fureby 
[ 6 ] points out that the numerical result is relatively in- 
sensitive to the LES subgrid model but more dependent 
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Figure 4: Instantaneous isobars at time step 740,000. 

on the grid sizes. Generally, a finer grid yields more de- 
tails in the acoustic field and higher SPL. It is important 
to choose a grid fine enough and an appropriate weighted 
average index number a. 

In the following, numerical results corresponding to 
flow speeds of 50 m/s and 26.8 m/s are presented and 
compared to the respective experimental data [15]. For 
the case with flow speed of 50 m/s, several consistency 
tests for the numerical work are also conducted. All the 
SPLs (sound pressure levels) are presented in dB based 
on their magnitudes against the standard 20 fiPa. 

4.1 Noise at flow speed of 50 m/s 

In order to have a sizable sampling time series for FFT 
(fast Fourier transform) analysis, with non-dimensional 
time step size At = 0.0025, over half million time steps 
are run. Figure 4 - Figure 5 demonstrate the isobars and 
v- velocity contours at time step 740,000. The unsteadi- 
ness and the large scale turbulence or instability waves 
are clearly displayed. In these figures, it is interesting to 
observe that the onset of the instability waves (vortices) 
does not occur right at the lip of the overhung tongue 
of the cavity, as expected for feedback oscillation, but 
occurs somewhat upstream. Evidently, the onset is trig- 
gered by the oscillation of cavity resonance. At the upper 
cavity surface around the overhung tongue, the turbulent 
boundary layer is disrupted and dominated by the large 
scale turbulence (vortices). 

With a laminar uniform ‘plug’ flow at the inflow 
boundary, due to the no-slip wall boundary condition, 
a boundary layer is naturally formed, developed, and 
turned turbulent near the overhung tongue of the cav- 
ity. No specified boundary layer thickness is given at 
the inlet. Table 1 demonstrates a comparison of the ex- 



Figure 5: Instantaneous v- velocity contours at time step 
740,000 . 


SPL in dB 



Figure 6: Sound pressure level (SPL) vs. frequency at a 
low speed of 50 m/s . 
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Table 1 : Comparison of experimental and computed tone noises at 50 m/s 


Exp. Freq. 
Hz 

Exp. SPL 
dB 

Comput. Freq. 
Hz 

Comput. SPL 
dB 

Compt. Freq. 
Hz 

lrg. domain 

Comput. SPL 
dB 

lrg. domain 

Noise origin 
? - unknown 

930 

103 

1075 

135 

1062 

131 

? 

1340 

107 

1326 

127 

1398 

130 

? 

1824 

134 

1837 

135 

1860 

129.6 

feedback 

2016 

113 

2062 

130 

2094 

119 

resonance 

2848 

106 

2962 

116 

2808 

111 

? 

3552 

111 

3513 

115 

3396 

106 

? 


perimental and computed data. 

From Table 1 , it is seen that for the fluid-dynamic feed- 
back tone, the numerical results agree very well with the 
experimental data, 1837 Hz. vs. 1824 Hz and 135 dB 
vs. 134 dB. in SPL, with a binwidth (or bandwidth) of 
47 Hz. The error in frequency or SPL is less than 1%, 
although Henderson stated that experimentally, ‘when 
multiple tones occur in the spectrum, the preferred or 
dominant mode often changes randomly. This can result 
in a change of 3 dB or more in the peak sound pressure 
level.’ [15]. For the Hemholtz cavity resonance, the fre- 
quencies are 2062 Hz vs. 2012 Hz with 2% error, while 
the error in SPL is larger possibly due to the reason just 
mentioned above. In addition, as shown in Table 1, mul- 
tiple tone peaks with lower SPL described by Henderson 
[15] as tones of unknown origins are also simulated. 

Some of these tone peaks of unknown origins in Fig. 6 
and Table 1 at the low frequency end seem to be subhar- 
monics, e.g. 1075 Hz vs. 2062 Hz. Subharmonics are 
often generated through vortex pairing, the excessively 
strong subharmonic spikes could be due to the 2-D sim- 
ulation of a 3-D flow. Such spikes are absent when a 
coarse grid is used (Fig. 7). 


SPL in dB, coarse grid 



Figure 7: Sound pressure level (SPL) vs. frequency at 
flow speed of 50 m/s with a coarse grid. 


sections. When comparing with the case of flow speed 
50 m/s in §4.1, it is referred to as the ‘standard case’. 


4.2 Consistency tests 

The same cavity noise problem is also computed with 
different grid size and larger computational domain di- 
mensions to ensure consistent and reasonable numeri- 
cal results are obtained. As mentioned above, we are 
trying to simulate the large scale turbulence itself. The 
more refined the grid, the more details of the turbulence 
emerge, as is experienced with a lot of LES (large eddy 
simulation) computations. In contrast to a deterministic 
flow computation problem, grid independence may not 
be tested in a rigorous way for the present noise problem. 
However, to a certain extent in a statistical way, some 
general agreements on frequencies may still be avail- 
able for different grids. Several other consistency tests 
such as running with a larger computational domain and 
longer time, as well as imposing a specified boundary 
layer thickness are also conducted in the following sub- 


4.2.1 Test of parallel computations The cav- 
ity noise problems are large scale computational prob- 
lems in the sense that a huge number of explicit time 
marching steps are needed. In order to reduce the turn- 
around time of the computation, parallel computation is 
adopted. The unstructured triangulated grid is first parti- 
tioned into subdomains by METIS [16], then MPI [17] is 
used to carry out data exchange between different CPU 
processors in the PC cluster. Sketch of the parallel com- 
putation can be found in [12, 13], 

Here, the same ‘standard’ cavity noise problem at flow 
speed of 50 m/s is tested on PC clusters with different 
number of processors ( 8 and 14) to confirm that the par- 
allel computations yield identical numerical results. 

4.2.2 Test with a coarse grid This cavity noise 
problem at flow speed 50 m/s was calculated with 
a coarser grid for grid independence test. Similarly, 
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SPL in dB, larger domain, 2,720,000 steps 


SPL in dB, 2.2cm B.L. thickness, 
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Figure 8: Sound pressure level (SPL) vs. frequency at 50 
m/s with larger domain and 2,720,000 marching steps. 


Figure 9: Sound pressure level (SPL) vs. frequency at 50 
m/s, with 2.2cm thick boundary layer imposed. 


nonuniform rectangular grid is used before it is triangu- 
lated. Although in the refined grid area Ax =0.00625 is 
about the same. Ay - .0114 is almost twice as large in 
the cavity opening area comparing to the current grid in 
use. Still, good agreement for the major modes with ex- 
perimental data and the current simulation was obtained. 
As shown in Fig. 7, the two shoulder points of the high- 
est spike correspond to 1875 Hz, 134.8 dB and 2080 Hz, 
137 dB. Since the binwidth (or bandwidth) is about 190 
Hz in this case, it is not clear that if they represent the 
fluid dynamic feed back oscillation or the fluid resonance 
mode, or both and form a broad peak. However, in terms 
of frequency, they are in good agreement with the ex- 
perimental data of 1824 Hz and 134 dB (fluid dynamic 
feedback) and the computed data of 1837 Hz and 135 dB 
with the ‘standard’ grid (Table 2). 

For the resonance mode, it is observed that the domi- 
nant frequency (2080 Hz) matches the experimental data 
well but the SPLs are higher (137 dB), the grid indepen- 
dence is thus somehow justified. But the tone peaks of 
unknown origins (broadband noises) are muffled, persis- 
tent higher harmonics are generated. 

4 . 2.3 Test on a larger computational domain 
and longer run The problem is also run with a large 
domain ( —2.8 ^ x ^ 4, 0 ^ y ^ 4). The grid size in 
this problem is about the same as in the ‘standard case’ 
in the critical area, i.e. the cavity interior and around 
the opening, but grows faster towards the inlet, outflow 
and top boundaries. Thus, a larger computational do- 
main is formed. With a time step size of At = 0.003, a 
much longer time of 2.72 million steps is run in attempt 
to find out the long term behavior of the simulation. The 


numerical results are presented in Fig. 8. Comparisons 
to the experimental data and the ‘standard case’ are still 
demonstrated in Table 1. Generally speaking, for the 
dominant oscillation modes, the aeroacoustic feedback 
and flow resonance, the error is about 2% which is quite 
acceptable, despite their lower SPLs as expected due to 
some coarse grid size. The performance is quite sim- 
ilar to the ‘standard case’ in that the subharmonic-like 
high SPL tones appear at the low frequency portion of 
the spectrum, possibly due to the 2-D simulation for 3-D 
flow. 

It is now demonstrated that, at least for the present 
simulation, the size of the computational domain (i.e., 
locations of the boundaries) and the long time running 
have little influence on the major numerical results. 

4 . 2.4 Tests with artificially imposed bound- 
ary layer thickness In all the above numerical simu- 
lations, a simple uniform plug flow is imposed at the inlet 
boundary. In this subsection, prescribed u-velocity pro- 
files are imposed at the inlet boundary to mimic the ac- 
tual turbulent boundary layers. Given a boundary thick- 
ness 6, the formal u-velocity profile at the inlet boundary 
is specified as: 

u = U 0 (y/S) 1/7 . 

In the following two test cases, the same computational 
domain and grid as in the ‘standard case’ are employed. 
In the first case, a thick turbulent boundary layer of 2.2 
cm is imposed at the inlet. After 1 .4 million time steps, 
the spectrum at the assigned point on the cavity vertical 
wall is depicted in Fig. 9. The prominent peaks corre- 
spond to 1660 Hz, 127 dB and 3320 Hz, 125 dB. Evi- 
dently, the second peak is the first harmonic of the first 


N AS A/CR— 2004-212892 


7 





Table 2: comparison of experimental and computed tone noises at 50 m/s, with coarse grid 


Exp. Freq. 
Hz 

Exp. SPL 
dB 

Comput. Freq. 
Hz 

Comput. SPL 
dB 

Compt. Freq. 
Hz 

coarse grid 

Comput. SPL 
dB 

coarse grid 

Noise origin 
? - unknown 

930 

103 

1075 

135 

- 

- 

? 

1340 

107 

1326 

127 

- 

- 

? 

1824 

134 

1837 

135 

1875 

134.8 

feedback 

2016 

113 

2062 

130 

2080 

137 

resonance 

2848 

106 

2962 

116 

- 

- 

? 

3552 

111 

3513 

115 

3960 

98 

? 



Figure 10: Sound pressure level (SPL) vs. frequency at Figure 11: Sound pressure level fSPL) vs. frequency at a 

50 m/s, with 1 .2 cm thick artificial boundary layer at inlet low speed of 26.8 m/s . 


peak. Note that the spectrum exhibits ample high fre- 
quency harmonics of a single tone, this implies numeri- 
cally, there exists a highly nonlinear wave. The pattern is 
quite different from the experimental one [15]. 

In the second case, the imposed boundary layer thick- 
ness at the inlet is reduced to 1 .2 cm. Figure 10 illustrates 
the spectrum after 1.4 million time steps. A slightly 
larger At = 0.003 is used. This time, the prominent 
peaks correspond to 1824 Hz, 121 dB and 3642 Hz, 124 
dB. It is also noted that in Fig. 10, to the left of the peak 
of 1824 Hz, there is a smaller peak at about 2090 Hz, 
which may correspond to the flow resonance frequency 
of 2016 Hz. In terms of frequency, the results are ex- 
ceedingly good, almost agreeing exactly with the exper- 
imental data for ‘thin’ boundary layer [15]. There is no 
surprise since <5 = 1 .2cm is exactly the boundary layer 
thickness from the experiment. However, in terms of 
SPL, the numerical ones are still low probably because 
the grid is still not fine enough. In addition, high SPL 


unwanted peaks appear at the low frequency end of the 
spectrum possibly due to 2-D simulation. 

The numerical experiments show that the prescribed 
boundary layer thickness has considerable effect on the 
spectrum and SPL. Without any knowledge or hint about 
the boundary layer, an arbitrarily prescribed boundary 
layer thickness may not work well and it is probably best 
to let the flow develop its own. 

4.3 Noise at at a lower flow speed of 26.8 m/s 

Numerical experiments were also conducted for an even 
lower flow speed of 26.8 m/s for the same cavity geom- 
etry. Table 3 lists a comparison of the experimental and 
computed data. As in the ‘standard case’, uniform plug 
flow is imposed at the inlet boundary. 

As in the previous case, from Table 3, it is observed 
that for the fluid-dynamic feedback tone, the numerical 
results agree very well with the experimental data, 1 146 
Hz. vs. 1168 Hz and 103 dB vs. 99 dB in SPL, with 
a binwidth (or bandwidth) of 95 Hz. The error in fre- 
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Table 3: Comparison of experimental and computed tone noises at 26.8 m/s 


Exp. Freq. 
Hz 

Exp. SPL 
dB 

Comput. Freq. 
Hz 

Comput. SPL 
dB 

Noise origin 
? - unknown 

928 

97 

856 

105 

? 

1168 

99 

1146 

103 

feedback 

1380 

95 

1420 

105 

? 

1890 

103 

1896 

105 

? 

1984 

101 

1988 

106 

resonance 


quency or SPL is about 3 — 4%. For the cavity resonance, 
the frequencies are 1988 Hz vs. 1984 Hz with 0.2% error, 
and the SPLs are respectively 106 dB and 101 dB. There 
are also other noise tones of unknown origin, the errors in 
frequencies and SPLs with these tones are larger. Sim- 
ilar to the case at flow speed of 50 m/s, their SPL are 
relatively high, possibly due to the 2-D simulation of a 
3-D flow. From Fig. 11, there is probably another tone 
of unknown origin - the one at 1896 Hz. Its frequency 
is so close to the resonance frequency 1988 Hz and the 
binwidth is 95 Hz, making it hard to distinguish the two 
frequencies. The two frequency spikes coalesce into a 
single broad peak. According to Henderson’s observa- 
tion from the experiments [15], well defined discrete fre- 
quencies do not appear consistently until the flow veloc- 
ity reaches 35 m/s. Based on this argument, the numer- 
ical results for the lower flow speed are indeed in very 
good agreement with the experiments. 

4.3.1 Test with a coarse grid The same prob- 
lem was tested on a coarse grid [18] similar to the one in 
§4.2.2. Due to lack of resolution, the SPL is much lower 
and only one spike at 1916 Hz appears in the spectrum, 
which may correspond to either the resonance frequency 
at 1984 Hz with 3.5%) error or the frequency spike of 
unknown origin at 1890 Hz. In terms of frequency, the 
approximation is acceptable. 

5 Concluding Remarks 

In this paper, a MILES type numerical approach is 
adopted to simulate the cavity tone noises arising from 
large scale turbulence. No explicit turbulence or subgrid 
models or acoustic models are used. As there are two 
dominant oscillation modes - aeroacoustic feedback and 
flow resonance competing with each other, accompanied 
by some tones of unknown origins, or broadband noise, 
it is important to select carefully the grid size and other 
parameters in the computations. Generally good results 
are obtained and compare favorably to the experimental 
ones, although 2-D simulation for 3-D flows may gen- 
erate some excessively strong subharmonics towards the 
lower end of the spectrum. 


In the MILES type approach, we attempt to simu- 
late only the large scale turbulence, which is believed 
to be the tone noise generating source. The smaller 
scale turbulences (higher frequencies and wave numbers) 
are practically ignored without an explicit subgrid scale 
model, since we are interested primarily in the aeroa- 
coustic data. The CE/SE method is a scheme for solving 
conservation laws. When the divergence theorem (2) is 
applied to find the cell average, it also plays a role as a 
filter ( such as the one used in LES approach). High fre- 
quency and high wave number disturbances are filtered 
out and aliasing errors are unlikely to occur. As pointed 
out by Fureby [6], choosing the parameters or limiters in 
the MILES approach, is somewhat equivalent to choos- 
ing a SGS model in the conventional LES approach. The 
choice is not straightforward in the presence of compli- 
cated flows. 
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